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Abstract. 



The competition between liquid and solid states of strongly correlated electron systems 



occurs in a novel way in a strong magnetic field. The fact that certain Landau level filling 
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factors are especially favorable for the formation of strongly correlated liquid states, gives 
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\0 ' rise to the fractional quantum Hall effect. In this article we discuss some consequences of 
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' for the thermodynamics of the liquid-solid transition. 
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the existence of incompressible states with fractionally charged quasiparticle excitations 
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g ■ 1. INTRODUCTION 



In the limit of strong magnetic fields electrons in two-dimensional (2D) systems are 
confined to the single-particle states of minimum quantized kinetic energy, i.e. to states 
in the lowest Landau level. In this limit all states of the system have the same kinetic 
energy so that the ground state, and for k^T <C hu c also all thermodynamic properties, 
are determined entirely by electron-electron interactions. In this article we restrict our 
attention completely to the strong magnetic field limit where only the lowest spin-polarized 
Landau level is relevant. 1 One single-particle state exists in the lowest Landau level for each 
magnetic flux quantum piercing the system and the electron density is usually measured in 
terms of the Landau level filling factor v = u^q/B = nhc/eB where n is the areal electron 



density and is the electron magnetic flux quantum. The lowest interaction energy state 
of 2D electrons is the Wigner 2 crystal state in which electrons are localized at the sites 
of a triangular lattice. At zero field it is the kinetic-energy cost of localizing electrons 
which prohibits the formation of the crystal state except at extremely low densities. In 
a strong magnetic field, liquid and crystalline states of electrons compete in a different 
way; electrons cannot be localized to a length smaller than the cyclotron orbit radius 
corresponding to the quantized lowest kinetic energy, £ = (hc/eB) 1 / 2 . Only when £ is 
much smaller than the typical distance between electrons, i.e. only when v <C 1 should the 
crystal state be expected to be the ground state. For weaker fields (larger filling factors) 
the ground state should be expected to be a liquid 3 ' 4 . 

The discovery of the fractional quantum Hall effect 5 '^ and its subsequent explanation 8 
established the intricate and exotic nature of the electron liquid states which occur in 
the strong magnetic field limit. More recently data from transport experiments 9 , and 
also from magneto-optical 10 and frequency-dependent conductivity experiments 11 , has 
accumulated 12 ' 13 which suggests the occurrence of the electron solid state. The exper- 
imental results are consistent with theoretical expectations that liquid states should tend 
to be favored near the filling factors where incompressible states occur and especially for 
v near 1/m where m is an odd integer. It is at these filling factors that the most stable 
incompressible liquid states, those first discovered by Laughlin 8 , occur. Many aspects of 
the experiments which seem to indicate a Wigner crystal state remain poorly understood, 
as discussed elsewhere in this volume, and it seems likely that in real samples disorder 
must play an important role 14 in determining the way in which Wigner crystallization is 
manifested. In this article we discuss only the ideal limit in which no disorder is present. 
We take the view 15 that in this limit the transition between fractional Hall liquid states 
and the Wigner crystal state is strongly first order so that we can study the thermody- 
namics of the phase transition by comparing separate estimates of the filling factor and 



temperature dependence of the free energies of liquid and solid states. In Section II we 
discuss the physics of fractional Hall liquid states for v near 1/m and estimates of the 
temperature and filling factor dependence of the free energy which result. In Section III 
we briefly review what is known about the filling factor and temperature dependence of the 
free energy of the Wigner crystal state. Our considerations are based primarily on the har- 
monic approximation for the lattice state. The filling factor and temperature dependences 
of the free energy are very different for the two states and as we discuss in Section IV 
these differences suggest the occurrence of several unusual features in the phase boudary. 
In Section V we briefly summarize our findings. 



2. THERMODYNAMICS OF THE LIQUID STATES 



For pairs of electrons restricted to the lowest Landau level the possible relative motion 
states may be labelled by angular momentum m = 0, 1, • • •. There is only one relative 
motion state of the pair for each value of m and for the fully spin-polarized states of the 
electron gas assumed here electrons can be found only in states of odd relative angular 
momentum. The Hamiltonian of the system can be completely specified in terms of the 
interaction energies of pair of electrons with relative angular momentum m, V m . For 
N$ = = m(N — 1) + 1 and any odd integer m it is possible to show 7 that there 

is only one state in the many-particle Hilbert space for which the amplitude for finding 
any pair of electrons in a state of relative angular momentum less than m is zero. (Here 
N is the number of electrons and A is the area of the 2D system.) That state is the 
strongly-correlated electron liquid state discovered by Laughlin, 8 ' 16 Since electrons 

are closer together when they are in a state of lower relative angular momentum it follows, 
for sufficiently short-range repulsive interactions, that at this particular magnetic field the 
Laughlin state will be separated from all other states in the Hilbert space by an energy 



gap A ~ V m -2- Similary the chemical potential at zero temperature will jump from a 
value ~ V m to a value ~ V m -2 when the filling factor increases beyond v = 1/m. The 
existence of the energy gap and the associated chemical potential discontinuity leads to 
the fractional quantum Hall effect. The chemical potential discontinuity must also lead to 
anomolies in the thermodynamic properties of the electron system at low temperatures for 
v near 1/m. 

For = m(N — 1) + 1 + the many-body states in which relative angular momenta 
less than m can be completely avoided can be mapped to the many-body states for N q h 
fermion holes in a Landau level with a degeneracy equal 17 to = N + — 1. The 
quasiholes must have fractional charge e/m since m of them are created when one electron 
is removed from the sytem at fixed magnetic fields. When the fractional quantum Hall 
effect occurs states in this class should be separated from higher energy states by a gap 
~ A. Similarly, it can be convincingly argued 18 ' 19 ' 20 from several different points of view 
that the low-energy states for = m(N — 1) + 1 — N qe can be mapped to the many-body 
states of N qe fractionally charged fermion particles in a Landau level with degeneracy = 
N + 1 — N qe . We refer to these excitations with fractional negative charge as quasielectrons 
and use the word quasiparticles to refer to the fractionally charged excitations generically. 
Note that when they are considered as fermions the quasiholes see both electrons and 
other quasiparticles as sources of effective magnetic flux with the flux from quasiparticles 
directed in opposition to that from the other sources. In the following we assume that this 
picture can be extended to apply to low-lying excitations as well as to the ground state. 
We estimate the thermodynamic properties of the electron system by assuming that its 
elementary excitations are quasiparticle-quasihole pairs. Since we are interested in filling 
factors near v = 1/m and low-temperatures we will assume that the quasiparticles are 
sufficiently dilute that we can neglect their mutual interactions. This approximation is 
adequate for most purposes 21 and as we shall see below some interesting non-trivial effects 



exist in the distribution functions for these ideal gases even when quasiparticle interactions 
are neglected. 

We consider the system at a constant magnetic field. Dropping terms ~ 1/N we define 
N = N — Nq where Nq = N^/m so that by total charge conservation N qe — N q h + mN. 
The energies of the states of the sytem at a given iV and hence at a given N may be 
expressed in terms of the number of quasielectrons and quasiholes: 

E(N qe , N qh ) = e iVo + e qe N qe + e qh N qh . (1) 

Here eg is the energy per electron of the incompressible state while e qe and are the 
energies to make quasiholes and quasiparticles at fixed magnetic field. For the physically 
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interesting case of Coulomb interactions between the electrons estimates exist for these 
three parameters: For m = 5 and in e 2 /£ units 23 €ql ~ —0.3277, « 0.1072, and 
e qe —0.076. (We use e 2 /£ units for energies throughout this paper.) Notice that the 
quasielectron-quasihole pair creation energy A = + e qp ~ 0.031. This means for = 
m(N — 1) + 1 the liquid ground state is separated from excited states in this approximation 
by an excitation gap equal 24 to A. At zero temperature the upward discontinuity in the 
chemical potential is mA at = m(N — 1) + 1. 

We evaluate thermodynamic properties at finite temperature in the grand canonical 
ensemble. The grand partition function is 

Z G = exp(-/3(e - fj)N ) ^ 5 N qa -N qht mN ex ^^ N ^l m ~ ^H^qe, N<f>, e qe )} 

N,N qe ,N qh 

e W [-^N qh /m - (3F(°\N qh , e qh )}, (2) 

where = N + N q ^ — N qe is the degeneracy of the quasiparticle and quasihole Landau 
levels. We see below that the dependence of the quasi-Landau-level degeneracies on the 
number of quasiholes and quasiparticles which are present on the chemical potential un- 
usual. In Eq. (2) F^(N, N^, e) is the finite temperature free energy for iV non-interacting 



Fermi particles in a system with a single energy level of energy e and degeneracy N^, which 
is known exactly 2 ^: 

F( \N, N+, e) = N^eu + k B T(uln(u) + ((1 - v) ln(l - v)\ = N^°\v, e) (3) 

where v = N/N^. To evaluate Eq. (2) we use the 5 function to eliminate the sum over 
N. At a given T and /U the most probable values of N qe and N q h can be determined by 
finding the extrema of the summand in Eq. (2). Setting the derivative with respect to N qe 
to zero gives 

H = m/jqe - (m - l)k B T[\n(l - v qe ) + ln(l - u qh )] (4) 

while setting the derivative with respect to N q ^ to zero gives 

/j = -mfi qh - (m - l)fc B Tpn(l - v qe ) + ln(l - v qh )]. (5) 

Here u qp = Nqp/N^ is the quasiparticle Landau level filling factor. The second terms on 
the right hand side of Eq. (4) and Eq. (5) come from the derivatives of the quasiparticle 
Landau level degenearacies with respect to the quasiparticle and quasihole numbers. In 
these equations we have defined 

Np = ^pEl^El = eqp + k B THv qp /{l - v qp )]. (6) 

so that the quasiparticle filling factors are related to their chemical potentials by the usual 
Fermi distribution function. Subtracting Eq. (4) and Eq. (5) we see that fx qe = —n q h so 
that , given fx, (i qe can be determined by requiring Eq. (4) to be satisfied when = —\x qe - 
The electron filling factor v = N/N^ can be related to the quasiparticle filling factors: 

m+ (m- l)(v qe - u qh )' 
We will be interested in comparing the Helmholtz free energy of the liquid at a particular 

filling filling factor with that of the solid. Noting that fluctuations become negligible in 

the thermodynamic limit we find from Eq. (2) that 

F = -k B T\nZ G + ^N = e iV + ^[/(°)(e ge , u qe ) + / (0) (e^, u qh )] (8) 



so that the free energy per electron is 

f=^T + - (rn - l))[f { °\e qe , v qe ) + fW(e qh , u qh )]. (9) 

1 1 Lis 



To determine v qe and v q ^ at a given filling factor and temperature we express the right 
hand side of Eq. (7) in terms of the quasiparticle energies and z qe = exp[/j, qe /(kBT)] = 
exp[— /^/(fc^T)] and solve for its value. Given z qe we can easily calculate v qe , v q ^ and 
the free energy. Some typical results are shown in Fig.[l]. The most noteworthy feature 
in this figure is related to the entropy, 

At v = 1/m the ground state is non-degenerate and separated from excited states by the 
gap A. It follows that the entropy vanishes like exp[— A/(2kgT)] at low temperatures 
so that the free energy, which must decrease monotonically with temperature, is nearly 
constant. For v ^ 1/m the ground state is degenerate as discussed near the beginning of 
this section, and the entropy approaches a constant. For v < 1/m, the degeneracy of the 
many-body ground state is 



N t N qh) 



N qh 



9 = 

so that the entropy at zero temperature in the thermodynamic limit is 
S 



= (1 + v~ l - m) ln[l + u- 1 - m] - (z/" 1 - m) ln^ 1 - m]. (12) 
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Similary the zero temperature entropy for v > 1/m is 



= (m — v ) ln[(m — v )' 



Nk B 

- (1 - 2(m - u- 1 )) ln[(l - 2(m - z/" 1 ))/(l - (m - (13) 

The finite zero-temperature free energy gives rise to the linear decrease in free energy with 
increasing temperatures which is evident at low temperatures in Fig.[l]. In a more realis- 
tic model the entropy would drop to zero on a temperature scale reflecting quasiparticle 



interactions , however these interactions would not greatly alter the comparisons between 
liquid and solid free energies which we make below. 



3. THERMODYNAMICS OF THE HARMONIC ELECTRON SOLID 

In this section we briefly review known results for the thermodynamics of the harmonic 2D 
Wigner crystal in the strong-magnetic-field limit. In using the harmonic approximation 
for the electron crystal we are ignoring the possibility of large departures of the electrons 
from their lattice sites. The use of this approximation is consistent with our assumption 
that the phase-transition is strongly first order since the harmonic approximation can then 
still be valid when the solid melts. We remark that the harmonic approximation is still 
useful for estimating thermodynamic properties at low temperatures even though true 
long-range-order in the Wigner crystal cannot exist at finite temperatures. 

Adding the Lorentz force to the classical equations of motion it is easy to demonstate 
that the harmonic phonon energies in the strong magnetic field limit are related to the 
zero-field energies 27 by 

e+(S) = hu c + J2 e l(Q)/ 2huj c (14) 
A 

and 

e-{q) = [\{ex{q)]/noo ( , (15) 
A 

(See for example Ref.[28]. Here is a zero-field phonon energy.) At long wavelenths 
the zero-field resonance frequencies are a purely longitudinal plasmon mode, e 2 ^(g) = 
27rft 2 ne 2 g/m* and a linearly dispersing transverse mode et = ctq. It follows that e+(q) = 
Uu c + e^ )l /2huj c while e-(q) = c t qe pl {q) /Tao c ~ q 3 / 2 . 

At general wavelengths the zero-field phonon energies can be expressed in the form 29 

aV^^ (16) 



where e(q) is a readily calculable pure number, q = n 1 / 2 ^, and n is the areal density of 
the electron system. It follows that at strong fields 

e+(q) = hu c + \{e 2 /^' 2 £ e~ 2 (q) = hu c + (e 2 /£)u 3 / 2 ~e + (q) (17) 

A 

and 

e_($) = (e 2 /^)^ 3/2 n e ~A(q) = (e 2 /V /2 Mq) (18) 
A 

Note that the low energy modes are independent of the electron mass and therefore must 
correspond to intra-Landau-level excitations of the system. For kgT <C %uj c the mean 
harmonic-osillator quantum number is zero for the high-frequency cyclotron mode so that 
the internal energy is given by E = Nhuc/2 + E mad + N(e 2 /£)v 3 / 2 e(t) where 

e W = ^5>+(q) + (MMo)/*) + l)e-(q)], (is) 

Q 

Emad 1S the classical point lattice energy, n{x) = (exp(x) + 1) _1 and t = k^T j {e 2 1 €)v^l 2 . 
Similarly the entropy has contributions only from the intra-Landau-level phonon modes 
and depends on temperature only through the dimensionless paramter t; S = Nkgs(t) 
where s(t) is readily evaluated numerically. In Fig. [2], Fig. [3], and Fig. [4] we plot results 
for the specific heat {Cy = Nk^e' (£)), the entropy and the Helmholtz free energy (F = 
Nhu c /2 + A^(e 2 /^)(-0.78213z/ 1 /2 + u 3 / 2 f(t)) respectively. Here fit) = e(t) - ts(t) and 
the second term in the Free energy is the Madulung energy. All quantities are evaluated 

an 

for the triangular electron lattice which has the lowest free energy over the temperature 
range of interest. In the next section we use this free energy to construct an estimate of 
the shape of the liquid-solid phase boundary near an incompressible filling factor. The 
free-energy of the solid, unlike that of the liquid, has a smooth dependence on filling 
factor. At low temperatures the vibrational contribution to the free-energy per electron 
is /(£) ~ — 9(e 2 /£)z/ 3 / 2 t 7 / 3 ; this power law follows from the </ 3 / 2 dispersion of the intra- 
Landau-level mode and should be compared to the linear temperature dependence which 



occurs except at v = 1/m for the liquid. Thus finite temperature effects should generally 
be expected to favor the liquid over the solid, as is commonly the case. An exception 
occurs for v very close to 1/m where there are few low energy excitations for the liquid. 

When leading anharmonic corrections are included the ground state energy of the 
Wigner crystal is given by 3 

e s = Nhu c /2 + N(e 2 /£) [-0.78213z/ 1/2 + 0.2410z/ 3 / 2 + 0.087z/ 5 / 2 ] (20) 

In Eq. (20) the first term in square brackets is the energy of the classical Wigner lattice, the 
second term is the harmonic zero point energy, and the last term comes from anharmonic 
corrections. (The second term should be compared with the zero-temperature limit of the 
vibrational free energy shown in Fig. [4]). We see that for v ~ 1/5, the region of filling 
factor we focus on below, anharmonic corrections are ~ 10~ 3 (e 2 /£) per particle. A similar 
uncertainty exists in estimates of the energy per particle of the incompressible state. Below 
we take the attitude that the difference between the ground state energies of liquid and 
solid states at v = 1/m cannot be determined with sufficient accuracy and may in practice 
be altered somewhat by changes in the quantum confinement width of the two-dimensional 
electron gas. We will treat this difference as a phenomenological parameter. 



4. THE LIQUID-SOLID PHASE BOUNDARY 



In Fig. [5] we plot the filling factor dependences of the free-energies for liquid and solid 
states obtained using the approximations discussed in the previous two sections. These 
curves were obtained assumming that the temperature dependence of the small anharmonic 
contribution to the energy of the solid can be neglected. The ground state energy of the 
liquid state has the cusp responsible for the fractional quantum Hall effect at v = 1/m. The 
ground state energy of the liquid is lower than that of the solid for a finite range of filling 



factors around v = 1/m. Well away from v = 1/m the degeneracy of the liquid ground 
state causes the free-energy to decrease more rapidly for the liquid than for the solid. We 
have constructed the phase diagram for our model of the liquid-solid phase transition by 
mapping out the filling factors at which the free-energy differences cross zero. The results 
are shown in Fig. [6]. There is qualitative agreement between these phase boundaries and 
those mapped out by various experimental probes 13 of the two-dimensional electron gas 
system in a strong magnetic field in this regime of filling factors. We emphasize that this 
level of agreement is achieved without any adjustable parameters. 

The temperature scale in Fig. [6] should be compared with the classical melting temper- 
ature of the Wigner solid which is T^f ~ 6x lCP 2 z/ 1//2 (e 2 /£/c£). Thus, in our approximation, 
the solid phase persists to temperatures well above the classical melting temperature. This 
behaviour seems to be improbable, although experimental fingerprints of the solid phase 
do also seem to persist beyond the classical melting temperature. The temperature scale in 
our figures is strongly sensitive to the quasiparticle-quasihole creation energy and could be 
adjusted downward by adjusting A. Our results are also strongly sensitive to the ground 
state energy difference between the liquid and solid states at v = 1/m. This difference 
is obtained by subtracting two quantities which have a relative difference of only several 
parts per thousand and must be regarded as highly uncertain. It may even be altered 
importantly by changes in the effective interaction between electrons due to differences 
in the electronic width of a hetero junction or quantum well between different samples. 
For that reason we plot in Fig.[7a-e] a series of phase diagrams obtained by arbitrarily 
adjusting the ground state energy of the liquid at v = 1/m so that is passes through our 
approximate value of the ground state energy of the solid at this filling factor. Members of 
this sequence of phase diagrams should also apply qualitatively, with suitable downward 
adjustments to the temperature scale, to the shape of the phase boundary for v near 1/m 
at larger values of m since the ground state energy difference of the two states is expected 



to cross zero 4 '*™ with increasing m. As the ground state energy of the liquid is increased 
the width of the filling factor interval where the liquid is stable at T = decreases. At the 
same time the T = entropy of the liquid at the position of the phase boundary decreases 
so that the effect of increasing temperature in expanding the width of the liquid stability 
interval is diminished. This is seen in Fig. [7a] and Fig. [7b] where the phase boundary lines 
become more vertical. When the ground state energy is raised still further the entropy of 
the solid along the phase boundary can increase above that of the liquid as temperature 
increases so that, in accordance with the Clausius-Clapeyron equation 31 , the width of the 
liquid stability interval decreases with temperature. As we see in Fig. [7c] and Fig. [7d] the 
stability interval can decrease to zero so that the solid becomes stable at v = 1/m over 
some finite temperature interval before melting again. Finally when the ground state at 
v = 1/m is raised above that of the the solid, as in Fig. [7e] , the solid will melt at a tem- 
perature which is minimized at v = 1/m. It seems likely that this is typically the situation 
for v = 1/9 and may often be the case for v = 1/7. The behavior at a particular filling 
factor in a particular system may depend on the degree of Landau level mixing i.e., on the 
electron density 30 , as well as on geometric details of particular systems which influence 
effective electron-electron interactions. 



5. SUMMARY 



In this article we have taken the view that the transition between Wigner crystal 
and fluid states in two-dimensions in a strong magnetic field is strongly first order. The 
thermodynamics of the phase-transition is then dominated by the same anomoly which is 
responsible for the fractional quantum Hall effect; namely the existence of incompressible 
states. When interactions between fractionally charged quasiparticles are neglected the 
entropy of the liquid has a finite contribution at zero-temperature which vanishes as the 



incompressible filling factor is approached. As a result the shape of the liquid-solid phase 
boundary changes qualitatively as the filling factors along the phase boundary approach 
the incompressible filling factor. When the ground state energies of the Wigner solid 
and the incompressible liquid are nearly identical the liquid can freeze and remelt with 
increasing temperature. 
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FIGURES 

Fig.[l] Free energy versus temperature in a fractional Hall liquid for u near 1/m with m = 5. 

Fig. [2] Temperature dependence of the specific heat of the harmonic Wigner lattice in the 
strong magnetic field limit. 

Fig. [3] Temperature dependence of the entropy of harmonic Wigner lattice in the strong mag- 
netic field limit. 

Fig. [4] Temperature dependence of the free energy of the harmonic Wigner lattice in the strong 
magnetic field limit. The dashed line shows the temperature dependence of the energy. 

Fig. [5] Free energy difference between liquid and solid states as a function of filling factor for 
a series of temperatures. The temper autres listed are in units of e^/ikg. 

Fig. [6] Phase boundary between Wigner solid and fractional Hall liquid states of the two- 
dimensional electron gas at strong magnetic fields near v = 1/5. 

Fig. [7] Phase boundaries between Wigner solid and fractional Hall liquid states of the two- 
dimensional electron for arbitrarily adjusted incompressible ground state energies. 



